Three-dimensional nodal superconducting gap in single crystals Ba(Fei_ x Ni x ) 2 As2 
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The London penetration depth, A, is directly related to the density, n s , of the Cooper pairs 
(A 2 oc l/n s ) and its variation with temperature provides valuable insight into the pairing mechanism. 
Here we study the evolution with doping of the temperature dependence of the in-plane (A a &) 
and out-of-plane (A c ) penetration depths in single crystals of electron-doped Ba(Fei_ cc NL c )2As2. 
As is the case for other pnictides, A(T) ~ T n over the whole doping range and this behavior 
extends down to at least T = T c /100, setting a very small upper limit on the gap minimum. 
Furthermore, in the overdoped regime: 1) the exponent n becomes substantially smaller than 2, 
which is incompatible with the models that explain power-law behavior to be due to scattering; 
2) the exponent n becomes anisotropic, with A C (T) showing a clear T- linear behavior over a large 
temperature interval. These findings suggest that in the overdoped regime the superconducting 
gap in iron-based pnictide superconductors develops nodal structure, which unlike in the cuprates, 
cannot be understood within a two-dimensional picture. 



Determining the symmetry of the superconducting gap 
function, A(k), is one of the milestones in understand- 
ing the pairing mechanisms of superconductors. In iron- 
based pnictide superconductors the symmetry of the gap 
is still a debated issue both theoretically and experimen- 
tally. Some early measurements, including angle - re- 
solved electron spectroscopy point contact 
spectroscopy [5j, penetration depth [6|, [3] and specific 
heat [8] were interpreted in a framework of fully gapped 
superconductivity. Others, such as NMR [9], penetration 
depth 0, El EE EDI and thermal conductivity 0, El 
have suggested a strongly anisotropic gap. More recent 
systematic measurements of thermal conductivity over 
the whole doping phase diagram in Ba(Fei_ :E Co ;E )2As2 
(FeCo-122) single crystals have shown a clear evolution 
towards a more anisotropic gap in the overdoped regime 
[l6[- Several other studies of the overdoped composi- 
tions have even suggested nodes in the superconduction 
gap H [TtJ • Other works have proposed a non- universal 
gap structure in different pnictide families [18]. Such a 
diversity of observations and interpretations is due to sev- 
eral factors, including a substantially three-dimensional 
character of the Fermi surface and the gap [19], competi- 
tion of superconductivity and other ordered states in the 
underdoped regime and a possible significant (pair break- 
ing) contribution from scattering. Additional complica- 
tions are due to complex chemistry and poor control of 
the composition and crystallinity, especially in RFeAsO 
(1111) and BaK-122 families of compounds. High qual- 
ity homogeneous single crystals are required to address 
questions related to the gap structure. 

In BaFe2As2, superconductivity can be induced by par- 
tial substitution of Fe with Co, Pd, Ni , R h and even 
combinations of Co and Cu 0, EH, Hi (IE Hi- Large 
homogeneous single crystals can be prepared over the 
whole phase diagram and the doping levels can be con- 
trolled and determined with high precision [2l], 0, 0] • 
In the overdoped regime, where superconductivity is not 



affected by the coexisting magnetic order (as is the case in 
the underdoped samples), the superconducting T c scales 
well with the number of dopant electrons, obtained by 
assuming a rigid band approximation [H Q . 

Among the experimental probes of superconductivity, 
the London penetration depth, A, is directly related to 
the density, n s , of the Cooper pairs (A 2 ex l/n s ), and 
its variation with temperature provides a valuable in- 
sight into the pairing mechanism. For a superconductor 
with a gap that is nearly uniform in momentum space, 
A(k) w Ao, the density of quasiparticles, n n , is exponen- 
tially small at low temperatures (n n = l—n 3 ) and A(T) is 
expected to saturate exponentially below T w T c /3. On 
the other hand, in cuprate high temperature supercon- 
ductors, the penetration depth shows linear temperature 
dependence (A ex T), consistent with a gap having line 
nodes [25[. 

In the present work, we have studied the evolution 
of both the in-plane (A a &(T)) and inter-plane (A C (T)) 
London penetration depths in single crystals of electron- 
doped Ba(Fei_ x NL E )2As 2 (FeNi-122). We have found 
that close to the optimal doping level, x op ^=0.046 [26[, 
\ a t)(T) is best described by the power-law, A a &(T) ~ 
T n , with n > 2, similar to previous results on FeCo- 
122 [12L [l3j]- In the overdoped regime however, the ex- 
ponent n decreases significantly below two and A a &(T) 
shows a much stronger temperature variation. An even 
more drastic evolution of A(T) is observed in the inter- 
plane penetration depth. On the under-doped side, A C (T) 
shows a tendency to saturation at low temperatures, 
while in the over-doped regime it is T— linear below 
T c /3. These results imply a doping-induced evolution 
towards a superconducting gap with three-dimensional, 
nodal structure in the overdoped regime. 

Single crystals of electron doped Ba(Fei_ £C NL E )2As2 
were selected from the growth described in Ref. 0. 
X-ray diffraction, resistivity, magnetization, magneto- 
optics and wavelength dispersive spectroscopy elemental 
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analysis (WDS) have all consistently shown good quality 
single crystals with a small variation of the dopant con- 
centration over the sample and sharp superconducting 
transitions [2^ ]. 

The penetration depth was measured using a tun- 
nel diode resonator (TDR) technique [13]. The sam- 
ple was situated inside an inductor with ac magnetic 
field H ac « 20mOe. Variation of A(T) change the sam- 
ple's screening ability, and thus the effective inductance 
and the resonant frequency. The frequency shift, A/, is 
proportional to AA. Details of the measurements and 
data analysis are described elsewhere [28]. For H ac || c, 
screening currents flow only in the afr-plane and Af is 
only related to the in-plane penetration depth, A\ ab , 
through A//A/ a6 = A\ ab /R ab , where Afg b is the to- 
tal frequency shift when sample is inserted into the res- 
onator and R ab is an effective sample dimension that ac- 
counts for the shape and the finite thickness of the sam- 
ple [28|, [29[ . When the magnetic field is applied along 
the afr-plane (H ac || afr), screening currents flow both in 
the plane and between the planes, along the c-axis. This 
situation is depicted in Fig. 03 In this case, Af 1 - contains 
contributions from both AA a 5 and AA C . For a rectangu- 
lar sample of thicknesses 2£, width 2w and length Af 1 - 
is given by Eq. (pQ) 

Af 1 - _ AA Qb AA C = AX mix 

Af 1 - t w ~ R b U 

where R b is the effective dimension [28| in this configura- 
tion and the excitation field H ac is parallel to the longest 
side / as illustrated on the top sketch in Fig. [3j Knowing 
AX ab from the measurements with H ac || c and the sam- 
ple dimensions, one can obtain AA C from Eq. (pQ). How- 
ever, because 2w > 4x 2t in most cases, Af 1 - is in general 
dominated by the contribution from AA a &. The subtrac- 
tion of AA C becomes therefore prone to large errors. The 
alternative, more accurate approach is to measure the 
sample twice [30]. After the first measurement with the 
field along the longest side / (H ac || Z), the sample is cut 
along this / direction in two halves, so that the width 
(originally 2w) is reduced to w. Since the thickness 2t 
remains the same, we can now use Eq. (pQ) to calculate 
AA C without knowing AA a 5. Note that the length / and 
width w are in the crystallographic ab— plane, whereas 
the thickness 2t is measured along the c— axis. In our 
experiments, both approaches to estimate AA C (T) pro- 
duced the same temperature dependence, but the for- 
mer technique had a larger data scatter as expected. We 
therefore only report AA C obtained by cutting and re- 
measuring the sample. 

All measured samples showed sharp superconducting 
transitions with AT < IK (determined as the onset of a 
diamagnetic signal), consistent with transport and mag- 
netization measurements [24]. Figure [TJ a) summarizes 
the T(x) phase diagram showing structural/magnetic 
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FIG. 1. (Color online) (a) Temperature-doping phase dia- 
gram of Ba(Fei_ x Ni x )2As2. The inset shows the power-law 
exponent n(x) obtained by fitting to X a b(T) = a + bT n from 
the base temperature up to T/T c =0.3. (b) AX a b(T) for differ- 
ent doping levels plotted vs. (T/T c ) 2 . The arrows mark the 
temperature below which n becomes less than 2. The curves 
are shifted vertically for clarity. 

and superconducting transitions. Figure 01b) shows 
the low-temperature (T < 0.3T C ) behavior of the in- 
plane penetration depth for several Ni concentrations. 
The data plotted versus (T/T c ) 2 are linear for under- 
doped compositions and show a clear deviation towards a 
smaller power-law exponent (below temperatures marked 
by arrows in Fig. [T](b)) for overdoped samples. While at 
moderate doping levels the results are fully consistent 
with our previous measurements in FeCo-122 [l2|, [l3j], 
the behavior in the overdoped samples is clearly less 
quadratic. It should be noticed that in order to suppress 
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FIG. 2. (Color online) Main panel: AX a b(T) for an overdoped 
sample with x = 0.072, T c = 7.5 K, (symbols) and the fit to a 
power- law a + 6T n , n — 1.62 (dashed line). Upper inset: Data 
taken with a 3 He cryostat and with a dilution refrigerator 
(DR) showing a good overlap. Lower inset: AX a b(T) on the 
normalized temperature scale comparing underdoped, x = 
0.033, T c = 15 K (UNDER), optimally doped, x = 0.046, 
T c = 19.4 K (OPT), and overdoped, x = 0.072, T c = 7.5 K 
(OVER), samples. 



T c by the same amount, one needs a two times lower Ni 
concentration compared to Co. In FeCo-122, the samples 
never reached highly overdoped compositions equivalent 
to x = 0.072 of Ni shown in Fig. [1(b). Therefore, Ni 
doping has the advantage of spanning the phase diagram 
with smaller concentrations of dopant ions, which may 
act as the scattering centers. The evolution of the ex- 
ponent n with x is summarized in the inset in Fig. [1(a). 
These values were determined by fitting the curves to 
AA a 5 = a + bT n over the whole temperature range shown 
in Fig. [1(b). If fitted over a temperature range with a 
smaller upper limit, the overdoped samples show an ex- 
ponent n significantly smaller than 2. 

For any superconducting gap structure, an exponential 
temperature dependence of A is expected roughly below 
Tmin ~ 0.3T c A min / A max . To estimate A min , we used a 
dilution refrigerator (DR), reaching a base temperature 
of about 60 mK. As shown in Fig. [21 in the overdoped 
sample the power-law behavior with n « 1.6 persists at 
least down to T/T c — 0.03, setting an upper limit for 
the in-plane A min w 0.1A max . The upper inset in Fig. [2] 
shows a good agreement between DR and 3 He measure- 
ments performed on the same sample. The lower inset 
compares samples of different doping levels. The slope 
(parameter b of the fit) correlates with the results found 
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FIG. 3. (Color online) Schematics of magnetic field penetra- 
tion in the case of H ac \\ I for the whole sample [A] and after 
cutting in half along length / [B]. (Main panel) The mixed 
penetration depth AXmixiT) before [A] and after [B] cutting 
for the overdoped sample x=0.072, T c = 7.5 K. Inset shows 
similar data for the underdoped x=0.033, T c = 15 K. 



in thermal conductivity measurements [16[. 

Until now, most studies have focused on the in-plane 
penetration depth. However, given the three-dimensional 
character of the band structure, it is imperative to sys- 
tematically study the inter-plane penetration depth as 
well. Figure [3] shows the effective penetration depth, 
^mix (see Eq. [T]), for overdoped, £=0.072 (main panel), 
and underdoped, x=0.033 (inset), samples before (A) and 
after (B) cutting in half along the longest side (/-side) as 
illustrated schematically at the top of the figure. Al- 
ready in the raw data, it is apparent that the overdoped 
sample exhibits a much smaller exponent n compared to 
the in-plane penetration depth, while underdoped sam- 
ples show a tendency to saturate below 0.13T C . Using 
Eq.[T]we can now extract the true temperature dependent 
AA C (T). The result is shown in Fig. [H for two different 
overdoped samples of the same composition, x = 0.072 
having T c = 7.5 K and T c = 6.5 K, and for an under- 
doped sample with x = 0.033 having T c = 15 K. Since 
the thickness of the sample is smaller than its width, we 
estimate the resolution of this procedure for AA C to be 
about 10 nm, which is much lower than 0.2 nm for AA a 5. 
Nevertheless, the difference between the samples is obvi- 
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less regime [35j. However, neither of these models can ex- 
plain a close to linear variation of the penetration depth 
and the anisotropy of the power-law exponent, n (unless 
scattering is significantly anisotropic). While scattering 
may still play an important role in determining the ex- 
ponent n, our results suggest a highly anisotropic nodal 
three-dimensional superconducting gap in the overdoped 
regime. In other words, not only is the gap not universal 
across different pnictide families [Hj], it is not universal 
even withing the same family over an extended doping 
range. 
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FIG. 4. (Color online) AA C (T) for the underdoped, x = 0.033, 
T c = 15 K (UNDER), and for two overdoped, x = 0.072, 
T c = 7.5 K and T c = 6.5 K (OVER), samples. Dashed lines 
are linear fits. 

ous. The overdoped samples show a clear linear temper- 
ature variation up to T c /3, strongly suggesting nodes in 
the superconducting gap. The average (between the two 
samples) slope is large, about d\ c /dT w 300nm/K in- 
dicating a significant amount of thermally excited quasi- 
particles. By contrast, in the underdoped sample the 
inter-plane penetration depth saturates indicating a fully 
gapped state. If fitted to the power-law the exponent in 
the underdoped sample 2 < n < 3, depending on the 
fitting range. 

At face value, our data suggest the development of 
a three-dimensional nodal structure in the overdoped 
regime in FeNi-122 crystals. Indeed, the penetration 
depth depends on the angular average of the response 
function over the Fermi surface. Nodes, if present some- 
where on the Fermi surface, will affect the temperature 
dependence of both components of A(T). However, the 
major contribution still comes from the direction of the 
supercurrent flow, thus placing the nodes in the present 
case at or close to the poles of the Fermi surface. The 
nodal topologies that are consistent with our experimen- 
tal results are latitudinal circular line nodes located at 
the finite k z wave vector or a point (or extended area) 
polar node with a nonlinear (A(#) ~ 6 P , p > 1) variation 
of the superconducting gap with the polar angle, 0. 

Another conclusion from the present study can be 
made with regard to the importance of scattering. A 
close to T 2 power-law behavior could be explained to 
be due to scattering in a nodal-gap scenario 
the proposed extended s± symmetry [33] with a sign 
change between different sheets of the Fermi surface [34| 
or strong pair-breaking that drives the system into a gap- 
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